Flat-band ferromagnetism induced by off-site repulsions 
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Density matrix renormalization group method is used to analyze how the nearest-neighbor repul- 
sion V added to the Hubbard model on ID triangular lattice and a railway trestle (t-t') model will 
affect the electron-correlation dominated ferromagnetism arising from the interference (frustration) . 
Obtained phase diagram shows that there is a region in smaller-f' side where the critical on-site 
repulsion above which the system becomes ferromagnetic is reduced when the off-site repulsion is 
introduced. 

PACS numbers: 71.10.-w,71.10.Fd,75.10.Lp 



The itinerant-electron ferromagnetism arising from 
electron correlations has a long history, dating back to 
the works by Hubbarctl, Kanamoril and Gutzwillera for 
the Hubbard model, ft has become increasingly clear 
that the criterion for the ferromagnetism is a formidable 
question when the correlation effect is fully taken into ac- 
count. There is a rigorous work by Nagaokau in the limit 
of infinite repulsion and infinitesimal doping, but the fer- 
romagnetism there is singular in Jjhat the spin stiffness 
vanishes as an inverse system sizeQ. 

Recently™ a new light has been shed on. the problem, 
when LiebQ, and later Mielke and TasakiD, have shown 
that flat (dispersionless) bands in the one-electron band 
structure are good news for the ferromagnetism. Re- 
markably, we can show that the correlation rigorously 
guarantees the ferromagnetism for arbitrary strength of 
the repulsion U when the flat band is half-filled. 

Occurrence of flat bands requires two classes of special 
lattice structures. Lieb's model exploits the lattices that 
have different numbers of sublattice sites, while Mielke's 
and Tasaki's model provide the flat bands from interfer- 
ences between the nearest-neighbor transfer t and more 
distant transfers t', which are assumed, to have sizeable 
magnitude (t! ~ t). Kusakabe and AokH have shown that 
the ferromagnetism is indeed stable, since, first, the spin 
stiffness is finite both in U/t — > and — > oo limits despite 
an apparent lack of relevant energy scales, and second-, 
the magnetism survives for finite dispersions. Tasakicl 
has shown exactly that the magnetism is not destroyed 
for some class of dispersive bands for sufficiently large U. 

These systems should be regarded as insulators be- 
cause half-filled flat bands are considered there. Whether 
a ferromagnetic ground state persists in metallic phases, 
i.e., for non-half-filled bands, is of great interest, since 
this will imply an itinerant ferromagnetism. Penc et aS3 
studied one dimensional(lD) models such as triangles 
connected linearly, which may be thought of as a strip 



cut out. from Kagome lattice, a realization of Mielke's 
modellld. Intuitively, a single triangle alone has a frus- 
tration in spin configurations, which results in an ef- 
fectively ferromagnetic exchange interaction. Penc et 
al have shown, using an effective Hamiltonian in some 
limit (U/t — > oo, etc), that the ID triangular lattice 
(Fig. lb) also has a ferromagnetic effective coupling be- 
tween spins. The system is then expected to show fer- 
romagnetism for general electron density in this limit. 
Sakamoto and Kubotj studied jtath the density matrix 
renormalization group (DMRG)El whether the ferr omag- 
netic ground state survives finite dispersions when hole 
is doped to the nearly flat band. 

On the other hand, Daul and NoackE-3 numerically de- 
termined the phase diagram of a ID Hubbard model hav- 
ing both nearest- and next-nearest-neighbor transfers (t- 
t' Hubbard model), which is topologically equivalent to 
a railway trestle (Fig.lc). They have found a large fer- 
romagnetic region for finite densities and finite on-site 
interactions. While a ID triangular lattice is a two-band 
system with the lower band being flat, the trestle has a 
single band, but its bottom can be nearly flat depend- 
ing on the value of t' . Hence we may include this ferro- 
magnetipn in the flat-band ferromagnetism in a broader 
context 

In the studies mentioned above, only the on-site 
repulsion(/7) is considered. In real materials, however, 
there are additional terms of appreciable strengths such 
as nearest- neighbor charge-charge interaction(y), bond- 
charge interaction(AT), exchange interaction(F), or ao.- 
site pair-hopping(_F'), which have been investigated!!^. 
For itinerant ferromagnetism, Kollar, Strack and Voll- 
hardt have derived sufficient conditions for realizing a 
ferromagnetic ground state for general lattices with one 
hole in a half-filled band (generalized Nagaoka's case). 
They conclude that F is important in stabilizing ferro- 
magnetism for finite U, while for the special case of X = t 
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the ferromagnetism is stable for F = 0. 

The purpose of the present paper is to study the ef- 
fect of the nearest-neighbor repulsion V for the flat-band 
ferromagnetism as exemplified by the ID triangular lat- 
tice or trestle. The reason why we focus on the effect 
of V is that V usually takes the largest value among 
V,X,F,F', playing an important role in such materials 
as organic compounds. Physically, the repulsion V, de- 
spite being a charge-charge interaction, can affect the 
magnetism through Pauli's exclusion that makes the in- 
teraction effectively spin-dependent, but how the effect is 
exerted is a non-trivial question. Here we have employed 
DMRG, a powerful method for investigating strongly cor- 
related ID systems, to obtain the phase diagram against 
U and t'. We have found that the ferromagnetic region 
shifts to the smaller-i' side when V is switched on. 

It is heuristic to start with a single triangle consisting 
of sites 1,2,3 (Fig{|a). The Hamiltonian is given by 
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in standard notations, where t is the transfer between 1,2 
and 2,3, t' the transfer between 1,3 with j = 1, and U 
the on site repulsion. The nearest-neighbor repulsion V 
is assumed to act between 1,2 and 2,3 only. Hereafter we 
take t = 1 as a unit of energy and consider the case of 
t' > 0, which favors the occurrence of flat bands. 

In FigJ|, we show for two electrons on a triangle the dif- 
ference in energy, A = E t — E s , between the lowest spin- 
triplet state (E t ) and the lowest spin-singlet state(-E s ) as 
a function of V for U — 10 with t' varied from 0.15 to 
0.25. Unexpectedly, the curves can be non-monotonic: 
For t' — 0.2 the ground state, which is a spin-singlet at 
V = 0, becomes a triplet around V = 5, then re-enters 
into a singlet for V > 7. 

The curious behavior can be understood intuitively as 
follows. Let us assume that the effect of the repulsion V 
can be taken into account by reducing t to a smaller t e g, 
because V reduces the amplitude of an electron residing 
on the top site 1. This is reminiscent of the approxima- 
tion often adopted in the t-J model, where the effect of 
the infinite on-site repulsion is taken into account by a 
reduction in t. When we vary t' in the absence of V, the 
high-spin state becomes most stable (A takes its mini- 
mum) at t' = t as shown in Fig.|| for U = 10. In the 
above argument V — 5 should then correspond to the 
case where t is reduced to t cS ~ t'(= 0.2 here). 

Now, an obvious question is: can this ferromagnetism 
induced by a repulsion survive when we connect the tri- 
angles into a ID chain? So we have applied the DMRG 
method to the Hubbard model on a ID triangular chain 
and to the trestle to look into such a possibility. 



The phase boundary between the ferromagnetic and 
paramagnetic phases can be determined in finite size 
systems as followsJi-3 For V ^ 0, we first calculate the 
ground-state energy (Epp) for fully spin-polarized state 
with DMRG. There the Hilbert space is much smaller 
than that of spinful electrons and the energy can be ob- 
tained very accurately. For V = 0, the ground state en- 
ergy can be readily obtained, since spinless fermions do 
not feel U. Next we calculate with DMRG the ground- 
state energy (-Eg) of a system having equal numbers of 
up-spins and down-spins (with S z — 0) to compare with 
E-pp. Alternatively we can directly calculate the total 
spin from 

Stot = 5^ (St ■ Sj). 

From these independent methods we can determine 
whether the ground state is ferromagnetic. In our cal- 
culation the two methods gave the same result. 

The calculation has been performed mainly for the 
number of atoms L = 39 with an open boundary condi- 
tion, and the convergence with respect to the system size 
has been confirmed by extending L to-59. We have kept 
up to 120 states per block at each stepE3 Using the finite- 
size algorithm, we swept the system about ten times. We 
stored the density matrix at each step to construct,good 
initial vector for each superblock diagonalizationLl We 
have obtained accurate wave functions with inverse it- 
eration and conjugate-gradient optimization after each 
diagonalization. 

Let us first look at the ID triangular chain, for which 
the one-electron energy band has a perfectly flat branch 
when t' = t/\/2. The Hamiltonian is given as before, 
where we have now 1 < i < L — 1, 1 < j < (L — l)/2, 
1 < k < L. We have found that the value of U required 
for the ferromagnetic ground state can indeed be re- 
duced for nonzero V for appropriate values of parameters, 
i.e., the band filling n=(number of electrons) /(number of 
sites) < 0.4 and t' < 0.5. A typical result is displayed in 
Fig.H, where we plot A = Epp — Eq as a function of U for 
n = 0.2 and t' — 0.31. The one-electron band structures 
are depicted in the inset of the figure. We can see that 
the ground state becomes ferromagnetic for U > 5 for 
V = 0, while the onset of the ferromagnetism is reduced 
to U = 4 for V = 2. Thus there is a parameter region in 
which the ferromagnetism is realized only when the off- 
site repulsion exists, as in the case of a single triangle. 

The result that this phenomenon occurs only for 
smaller band fillings (n < 0.5) may be understood as 
follows. In the above argument for a single triangle, we 
have assumed two electrons. If we envisage that we re- 
quire at least one hole per triangle (i.e., one hole per unit 
cell that consists of two sites) for the triangle lattice in a 
similar manner, this precisely amounts to n < 0.5. 

Let us next consider the t-t' Hubbard model(fig. |I|c). 
The Hamiltonian is given as 
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nim + i. 



In Fig.|| we again plot A = Epp — Eq as a function of U 
for t' = 0.16, n — 0.4, and V — or 4. It can be seen 
that the critical value, U Cl above which the ferromag- 
netism appears is reduced from U c = 7 down to U c = 5 
for V = 4. When we make the off-site repulsion too large 
(V = 20), the ferromagnetic ground state is washed away 
at least for U < 8. Thus, as in a triangle and in a triangle 
chain, an intermediate value of V makes the ground state 
to be ferromagnetic for appropriate values of t' , n and U. 

Figure || is the full phase diagram on the U — t' plane. 
From this we can see that the ferromagnetic region, which 
has a concave boundary for V — 0, shifts to the lower 
side of t' when V is introduced. The reduction of U c 
found above is one manifestation of this. Intuitively, the 
reason why we have a concave ferromagnetic boundary 
with a minimum around t' ~ 0.36i for V — is that the 
ferromagnetism becomes most favorable when the one- 
electron band, displayed in the accompanied panels in the 
figure, approaches a flat band. Namely, the band in the 
t — t' model cannot become exactly flat, but can become 
nearly so for t' ~ 0.36t, which is in between the single- 
minimum band and the double- minimum band. The shift 
of this concave curve with V may be understood qualita- 
tively from the above argument that V effectively reduces 
0.36i into 0.36t eS - 

For the generalized Nagaoka's ferromagnetism, a suffi- 
cient condition for a reduction of U c has been shown by 
Kollar et aS3, where for F ^ or X = t the condition 
on U is relaxed when we take an appropriate value of 
V. The present result implies that for small t 1 , off-site 
repulsion can stabilize the flat-band ferromagnetism as 
well. 

To summarize, we have found for (nearly) flat bands 
arising from frustration (ID t-t' model and triangular 
lattice) a region in the phase diagram where the criti- 
cal on-site repulsion for ferromagnetism is reduced when 
we introduce an off-site repulsion. Recently, ferromag- 
netism in three-dimensional-lattice structures (bcc, fee, 
etc) has received attentiorO ci. Since network formed 
by the transfers in bcc and fee lattices are also frustrated 
in that they comprise triangles, it may be interesting to 
analyze the effect of V in these systems. 

R.A. is grateful to K. Kusakabe and H. Sakamoto for 
useful discussions. Computations were done on SR2201 
at the Computer Center of the University of Tokyo and 
FACOM VPP 500/40 at the Supercomputer Center, In- 
stitute of Solid State Physics, University of Tokyo. 
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FIG. 1. A single triangle (a), ID triangular lattice(b), and 
the railway trestle (or t-t') model(c). 



FIG. 2. The difference in energy, E t — E s , between the 
lowest triplet state and the lowest singlet state for a single 
triangle plotted as a function V with U = 10 for t' varied 
from 0.15 to 0.25. 



FIG. 3. The difference in energy, Et — E s , between the 
lowest triplet state and the lowest singlet state for a single 
triangle plotted as a function t' for U = 10 with V = 0. 
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FIG. 4. The difference in energy (A) between the lowest 
ferromagnetic state and the ground-state for a ID triangular 
lattice plotted as a function of U for t' = 0.31, n = 0.2 with 
V — or 2. Solid circles represent the result for 39 sites, dia- 
monds for 59 sites, while the arrows indicate U c - Inset shows 
one-electron band structures, where the horizontal dashed line 
indicates the Fermi level. 



FIG. 5. A similar plot as in the previous figure for the t-t' 
Hubbard model for t' = 0.16, n = 0.4 with V = or 4. 



FIG. 6. The phase diagram (P: paramagnetic, F: ferro- 
magnetic) against U and t' for the t-t' Hubbard model for 
n — 0.4 with V = or 4. One-electron band structures are 
displayed in the accompanied panels for three typical values of 
t' labeled with (a),(b),(c), where the horizontal dashed lines 
indicate the Fermi level. 
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1D triangular lattice 



t'=0.31, n=0.2 
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t-V Hubbard model 
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